Computational scheme for ab-initio predictions of chemical compositions interfaces 

realized by deposition growth 
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We present a novel computational scheme to predict chemical compositions at interfaces as they 
emerge in a growth process. The scheme uses the Gibbs free energy of reaction associated with the 
formation of interfaces with a specific composition as predictor for their prevalence. It explicitly 
accounts for the growth conditions by rate-equation modeling of the deposition environment. The 
Bell-Evans-Polanyi principle motivates our emphasis on an effective nonequilibrium thermodynamic 
description inspired by chemical reaction theory. We illustrate the scheme by characterizing the 
interface between TiC and alumina. Equilibrium thermodynamics favors a nonbinding interface, 
being in conflict with the wear-resistant nature of TiC/alumina multilayer coatings. Our novel 
scheme predicts that deposition of a strongly adhering interface is favored under realistic conditions. 
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I. INTRODUCTION 

Density functional theory calculations are today rou- 
tinely applied to characterize structural and electronic 
properties of condensed matter systems. They serve as 
an important complement to and extension of experimen- 
tal methods [lj]. Atomic (and electronic) structure and 
chemical composition are inseparably interwoven. Infor- 
mation about chemical composition at, for example, in- 
terfaces (including surfaces) is therefore of great impor- 
tance for reliability of such calculations. 

Traditional ab initio atomistic thermodynamics meth- 
ods aim at describing and predicting compositions at ox- 
ide surfaces @, H| . These schemes assume that equilib- 
rium is established between the oxide and surrounding 
O2. This criterion is often justified for oxide surfaces 
that are in direct contact with an O-rich environment. 
For (solid-solid) interfaces, the situation is more compli- 
cated. Interfaces are typically exposed to the surround- 
ing environment only at the moment of creation. 1 Fur- 
thermore, oxides are seldomly grown directly from O2, 
and it is not clear what gas(es) they are in (dynamic) 
equilibrium with during deposition, if at all. A general- 
ization of the ab initio atomistic thermodynamics scheme 
from surfaces to interfaces [J] can therefore be problem- 
atic. In fact, we have shown Q that the equilibrium 
configuration at the interface between TiC and alumina 
predicted by such a generalized scheme does not describe 
the wear-resistance of TiC/alumina multilayer coatings 
We have attributed this inconsistency to the fact 
that the scheme does not account for the actual growth 
conditions. This may be a serious shortcoming also for 
other interfaces realized by deposition growth. 

This paper presents a computational scheme that ex- 
plicitly accounts for deposition conditions. At the same 
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1 We assume that the interface is created by a deposition process 
from a gas phase. 



time no a priori equilibrium assumptions are introduced. 
The scheme is therefore capable to predict chemical com- 
positions at interfaces (including surfaces) as they arise 
in a deposition process. While kinetic accounts are gener- 
ally required to describe the full nonequilibrium nature of 
growth, the Bell-Evans-Polanyi (BEP) principle Q mo- 
tivates an effective thermodynamic description inspired 
by chemical reaction theory. 

The key elements of our scheme are the use of Gibbs 
free energies of reaction as a predictor for the chemi- 
cal composition and a modeling of the deposition con- 
ditions in terms of rate equations describing the depo- 
sition environment. Our approach extends the method 
of ab initio thermodynamics of deposition growth for 
surfaces |8] to interface modeling. Applying the scheme 
to the TiC/alumina interface in a chemical vapor depo- 
sition (CVD) environment we predict the formation of 
strongly adhering structures in agreement with the wear- 
resistance [f| of CVD TiC/alumina multilayers. 

The paper is organized as follows. Section [TTI motivates 
the use of the Gibbs free energy of reaction as a predictor 
of chemical composition. Section UTTl contains the details 
of our modeling. Results are presented in Sec. IIVI and 
Sec. IVl contains our conclusions. 



II. PREDICTOR FOR AS-GROWN CHEMICAL 
COMPOSITION 

We use the Gibbs free energies of reaction G r as a pre- 
dictor for the prevalent chemical composition at surfaces 
and interfaces realized by deposition growth. We justify 
the use of G r as predictor using chemical reaction theory 
HE! and the BEP principle 0. 

A description of nonequilibrium growth generally re- 
quires a kinetic description. The prevalent surface or 
interface composition will in principle depend on the 
reaction barriers for all relevant processes. However, 
the BEP principle Q motivates our emphasis on an ef- 
fective nonequilibrium thermodynamic approach. The 
BEP principle states that the more exothermic a reac- 



tion is, the smaller is the associated reaction barrier. The 
nonequilibrium system will therefore continuously evolve 
along paths that release the largest Gibbs free energy of 
reaction. In some cases we can even obtain quantita- 
tive predictions for the prevalence of competing chemical 
compositions without having to calculate specific rates. 

The Gibbs free energy of reaction is defined as the 
gain in Gibbs free energy in a chemical reaction. For a 
general reaction, according to chemical reaction theory, 
G r is related to the forward (/) and backward (b) reaction 
rates Tf = kfH r [R r \ Vr and T b = k b Ii p \P p ] v p via 



-pG r = ini7/r b 



(1) 



Here f3 is the inverse temperature in units of energy, &{/,&} 
is forward and backward rate constant, [R r ] ([Pp]) is 
the concentration of the r-th (p-th) reactant R (product 
P), and ^{p, r } is the corresponding stoichiometric coeffi- 
cient. 2 The key observation is that the ratio of forward 
and backward rates in a given chemical process is en- 
tirely given by the nonequilibrium Gibbs free energy of 
reaction. This quantity can be evaluated without know- 
ing the kinetic details of the process. 

To illustrate the approach, we first discuss the prob- 
lem of surface termination in deposition growth [8| , seek- 
ing to predict which surface (A terminated or B termi- 
nated) will emerge when growing a binary material AB, 
see Fig. [TJ We consider two coupled reactions I and 77 
that change the surface termination from A to B and vice 
versa, 



I : B + React/ ^ A + Prod/ 
rf 



II : A + React// ^ B + Prod//. 
Tl 1 



(2a) 
(2b) 



Here, React{/ //} and Prod^/ //} collectively denote the 
reactants and products in the two reactions. The for- 
ward and backward rates are the same as in ([1) but have 
an additional index to differentiate between the two re- 
actions. Taking the difference between the free energies 
of reaction in I and 77, we find after exponentiation and 
use of the geometric mean {x,y) gm — (xy)^ 1 !" 1 
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The reactions in @ can also be described by the rate 
equations, 



d t P A = 
d t P B 
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FIG. 1: Illustration of two different terminations at the sur- 
face of a binary material AB. 



Here we have introduced the probability P{a,b} for ob- 
serving either an A or B terminated surface. Using the 
arithmetic mean (x,y) am — (x + y)/2, the steady-state 
solution provides a ratio between these probabilities, 



n 1 ) 



p B (n + ry) (ri,ry) a 



(5) 



In dynamic equilibrium [j| [13] , being characterized by 
Gl + G 1 / = or equivalents Tj/T^ 1 
by combining <j3j> and (0 that 



Y[/Tf, we find 



Pa> 



B 



exp 



(-P[Gl - ]/2) 



(G) 



In particular, if the reaction that creates an A-terminated 
surface has lower G r than the reaction that creates a 
B-terminated surface, here G> < G^ 1 , an A-terminated 
surface is more likely. Away from dynamic equilibrium, 
this relation is no longer exact. However, comparison of 
G r for different terminations still indicates which termi- 
nation is most prevalent, since the geometric and arith- 
metic means [Eqs. © and ©] are approximately equal, 
Ref. Q. 

By analogy to the case of surface terminations, we can 
also use differences in G r , see ©, as measure to predict 
chemical compositions at interfaces formed by deposition 
growth (and assumed to have retained the structure spec- 
ified by the deposition environment). Considering related 
interface configurations of similar thickness but different 
compositions, the one with lowest (most negative) Gibbs 
free energy of reaction is most likely to describe the na- 
ture of the interface as it is formed by deposition growth. 
We note that G r will in general decrease (if the reaction 
is favorable) as the number of constituents in the film 
describing the interface increases by any integer multiple 
of the bulk stoichiometry. The predictor is useful when 
comparing films of similar thickness. Again it is possible 
to obtain quantitative predictions without having to con- 
sider kinetic details of the deposition, see also Ref. [I lj . 



III. MODELING 



A. Materials Background 



2 In this section we count stoichiometric coefficients positively, op- 
posed to standard chemical nomenclature, where stoichiometric 
coefficients of reactants are counted negatively. 



We illustrate our computational scheme by studying 
the interface composition between TiC and alumina. 
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FIG. 2: Illustration of chemical vapor deposition of alumina. 
A H2-AICI3-CO2 gas mixture is supplied to a chamber at rate 
Rs- The chamber is kept at fixed temperature T and fixed 
total pressure p. The gases react to simultaneously form water 
and alumina. Reaction products and unused reactants are 
exhausted at rate Re- 




FIG. 3: Exfoliating (left panel) and wear-resistant alumina 
(right panel) films on TiC. Color coding: Ti = dark green, C 
= light green, Al = blue, O = red. 



TiC/alumina multilayers are commonly used as wear- 
resistant coating on cemented-carbide cutting tools @- 
They are routinely fabricated by chemical vapor deposi- 
tion (CVD). 

Figure [5] illustrates the experimental setup for CVD 
of alumina 12]. A H2-AICI3-CO2 supply gas mixture is 
injected into a hot chamber which is kept at a fixed tem- 
perature and a fixed total pressure. The CVD process 
proceeds in two steps which, however, take place in par- 
allel. Water forms at an (unknown) rate i?H 2 o according 
to 



->H 2 + CO. 



H 2 + C0 2 



Alumina is deposited at an (unknown) rate i?Ai 2 o 3 
cording to 



2AICI3 + 3H 2 



h A1 2 3 + 6HC1. 



(7) 



ac- 



(8) 



B. Ab initio structure search for thin alumina films 
and TiC/alumina interface models 

In Refs. 

nmm we have ( 

in collaboration with oth- 
ers) presented a computational strategy for identifying 
energetically favorable geometries of thin-film alumina 



on a TiC substrate. We have constructed alumina can- 
didate configurations from a pool of structural motifs as 
they exist in stable and metastable bulk alumina phases. 
The motifs are characterized by stacking of the O layers 
and by the coordination of the Al ions. We have consid- 
ered alumina films of thicknesses up to four O layers. In 
addition, we allowed for off-stoichiometries in these films. 

Each candidate configuration was structurally opti- 
mized by ab initio total-energy and force calculations 
and subsequent relaxation (using standard quasi-Newton 
and conjugate gradient algorithms). For a given thick- 
ness and stoichiometric composition, the optimal geom- 
etry was identified as the one with lowest total energy 
after structural relaxations. 

Figure [3] details the atomic structure of two optimized 
alumina films with a thickness of four O layers and dif- 
ferent stoichiometric compositions. These films arc re- 
ferred to as exfoliating (left) and wear-resistant (right) 
films. These terms refer to the adhesion properties of the 
two films (exfoliating = no or very weak adhesion, wear- 
resistant = strong adhesion). Details concerning adhe- 
sion properties of these films are given elsewhere pd| . 
For other thicknesses, the films possess similar structural 
(and adhesion) characteristics. Focusing on the films in 
Fig. [3] therefore does not pose a restriction of generality 
on the presented discussion of the TiC / alumina interface. 

The results of our (total-energy) structure search (for 
each class of AImOjv films) can be combined with the 
equilibrium thermodynamics of Refs. [2-4]). We find that 
this approach erroneously identifies the exfoliating film 
as thermodynamically stable This incorrect equilib- 
rium prediction applies to a wide range of temperatures 
and 2 pressures and is in appearand conflict with wear- 
resistance of TiC / alumina multilayer coatings @ . Below 
we show that a more physical account of alumina growth 
emerges with our new nonequilibrium thermodynamics 
theory. 



C. Formation of excess atoms and free energies of 
reaction 

The films shown in Fig. [3] can be divided into a stoi- 
chiometric part and an excess part. For a general AImOjv 
film, the number of stoichiometric units in the stoichio- 
metric part of the film is n a i U mina = min([M/2], [N/3]), 
where [x] is the largest integer smaller or equal to x. 
The corresponding number of excess Al and O atoms is 
An A i = N - 2n a i umina and An = M - 3n a i umina . 

We assume that the deposition of the stoichiometric 
parts of the films is described by (JSJ. Excess O can be 
deposited as 



C0 2 - 
H 2 



O c > 

► Cv 



CO, 
hH 2) 



(9a) 
(9b) 



depending on which reaction gives a lower free energy of 
reaction. (Excess Al is not considered here.) The free 
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energy of reaction is defined accordingly, 

qAI m o n _ G TiC / AlAf0jv - G T ic + "-M2O3 (6^hci - 2 maici 3 
+ An max[/i C0 -Mco 21 Mh 2 -Mh 2 o]- (10) 

D. Approximations for evaluation of G r 

We make standard approximations for the free energy 
of solid and gaseous constituents. The free energy of the 
TiC/alumina systems is approximated by their DFT to- 
tal energies, that is, G^oM ~ -E'soiid- 3 Vibrational effects 
are not considered [D, y§]. For gaseous constituents, we 
employ the ideal-gas approximation, 

IM(T,pi) = e t + A S °(T) + k B T\n{ Pi /p°). (11) 

Here 6j is the DFT total energy [l5[ of the gas phase 
species (molecule), and k B is the Boltzmann constant. 
A°(T) is the temperature dependence of fii at a fixed 
pressure p , and available for p = 1 atm in thermo- 
chemical tables fill ]. 

E. Rate-equation description of CVD of alumina 

The actual evaluation of the individual chemical po- 
tentials of the different gaseous constituents in the CVD 
chamber requires the specification of the associated par- 
tial pressures. In the actual fabrication process, these 
are, however, not directly controlled; only the tempera- 
ture and the total pressure are controlable. 

We describe the CVD process in terms of rate equa- 
tions for the individual (ideal gas) pressures, 

d tPl oc Ci Rs - -Re + vf 2 °Rn 2 o + ^ Al2 ° 3 i?Ai 2 o 3 - 

(12) 

Here, pi = Pi(t) is the pressure of chemical species i at 
time t inside the reaction chamber, p = p(t) = ^2 iP i(t) is 
the corresponding total pressure, a is the concentration 
of the chemical species i in the supply gas, and vf^ 2 ° 
and i/^ l2 ° 3 are the stoichiometric coefficients 4 of the the 
chemical species i in reaction and ©, respectively 
The rate at which the gas is supplied to (exhausted from) 
the chamber is Rs (Re) and the reaction rates for water 
production and alumina deposition are i?H 2 o an d -Rai 2 o 3 - 




FIG. 4: Gibbs free energies of reaction G r for exfoliating and 
wear-resistant alumina film (see Fig. [3]) as functions of the 
scaled reaction rate Th 2 o- The vertical line limits rn 2 o to 
the right and corresponds to dynamic equilibrium in J7|l. For 
simplicity, we assume J?ai 2 o 3 = -Rh 2 o/3. The difference in 
G r is a predictor for the relative presence of exfoliating and 
wear-resitant films, see ©. 

We use the resulting steady-state pressures 

1 H a O , A1 2 3 

_ Cj±rB 3 QVj +r A i 2 o 3 ^ , , 

Pi — P 1 , ) 

1 + ^A1 2 3 

as input for the evaluation of the individual chemical po- 
tentials (1TT1) . In (IT3)) we have introduced the scaled re- 
action rates r H2 o = Rr 2 o/Rs and r Al2 o 3 = Rai 2 oJRs- 

F. Limits on reaction rates 

The reaction rates i?H 2 o an d Ra\ 2 o 3 cannot assume 
any arbitrary value. For example, we require that 
-Rai 2 o 3 < -Rh 2 o/3 = Rm£) 3 , since reaction ^ requires 
three units of H2O, while reaction ([7]) produces only one. 
Additional constraints on this particular growth process 
and interface formation are discussed in Ref. 111. 



3 In fact, we correct the total energies of the films by subtract- 
ing the strain energy of the stoichiometric part of the film, 
-Efilm -» Eftlm ~ «A1 2 3 ^strain, where A strain is the is differ- 
ence between the strained (due to the expansion to the TiC lat- 
tice in the interface plane) and the unstrained bulk alumina per 
stoichiometric unit. 

4 Here we use the standard conventions that stoichiometric coeffi- 
cients are counted negative if a species is consumed and positive 
if a species is produced in a reaction. 



IV. RESULTS 

In Fig.|3]we plot the free energies of reaction for the two 
films shown in Fig. [3] as functions of the scaled reaction 
rate th 2 o- Deposition parameters (supply gas composi- 
tion, total pressure, and deposition temperature) as spec- 
ified in Ref. [HJ have been used. In this illustration of 
the approach we assume for simplicity i?Ai 2 o 3 = ^h 2 o/3- 
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Furthermore, the scaled reaction rate rn 2 o is limited to 
the right by the limit of dynamic equilibrium in (|7|). in- 
dicated by the vertical line. 

We find that, over the entire range of the reaction rate 
r H 2 Oi it is more favorable to grow wear-resistant overlay- 
ers than exfoliating overlayers. 5 Qualitatively, this result 
is independent of the choice of i?Ai 2 o 3 - We emphasize 
that this prediction is consistent with industrial use of 
TiC /alumina as wear-resistant coating; predictions based 
on an adaption of equilibrium-thermodynamics analysis 
0-0] to this system [j| are not. 

V. SUMMARY & CONCLUSIONS 

We have presented a novel computational scheme 
to predict chemical compositions at interfaces as they 
emerge in a growth process. The scheme uses the Gibbs 
free energy of reaction associated with the formation of 
interfaces with a specific composition as predictor for 
their prevalence. We explicitly account for the growth 
conditions by rate-equation modeling of the deposition 
environment. 

An earlier study of TiC/alumina interfaces [5] doc- 
umented that the predicted composition at this in- 
terface is in conflict with the wear-resistance of mul- 
tilayered TiC/alumina when using an equilibrium- 
thermodynamics scheme 041] ■ Our results demonstrate 
that a careful account of deposition conditions in inter- 

5 We notice that there is a third stoichiometric composition (not 
discussed in this work) that may be stabilized if t"h 2 ' s suffi- 
ciently far away from the dynamic equilibrium limit. The corre- 



face modeling is crucial for understanding the adhesion 
at CVD TiC/alumina. 

We expect that a similar analysis will be necessary also 
for other buried interfaces that form during a deposi- 
tion process in an environment that strongly differs from 
ambient conditions. We emphasize the predictive power 
of the here-presented method, adding to the DFT-based 
toolbox for accelerating innovation [17] . In principle, it 
allows for the determination of deposition conditions re- 
quired to experimentally create interfaces with a prede- 
termined compositions. We argue that combining the 
structure search method of Ref. [5[ and the new ab initio 
thermodynamics with, for example, a genetic algorithm 
18, 19] presents a valuable tool for characterization of 
surface and interface growth. 

Acknowledgments 

We thank Carlo Ruberto and Gerald D. Mahan for 
valuable discussions. We gratefully acknowledge support 
from the Swedish National Graduate School in Materials 
Science (NFSM), from the Swedish Foundation for 
Strategic Research (SSF) through ATOMICS, from the 
Swedish Governmental Agency for Innovation Systems 
(VINNOVA), from the Swedish Research Council (VR) 
and from the Swedish National Infrastructure for Com- 
puting (SNIC). 

sponding structure has also wear-resistant properties [Tl| . 
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